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An analytic model for non-spherical lenses in covariant Modified 
Newtonian Dynamics 
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1 INTRODUCTION 



The M odified Newtonian Dynamics paradigm (MOND, iMilgroml 
Il983h has been originally proposed as an alternative to galactic 
cold dark matter (CDM). This paradigm postulates a modification 
of Newton's laws in order to explain the rotation curves of spiral 
galaxies, by saying that the MONDian gravitational accelerations 
^M <K Qo ~ lO^'^ms"^ approach (gA-ao)''^ where gjv is the usual 
Newtonian gravitational acceleration. Nowadays, this paradigm is 
still better suited than dark matter to explain the observed conspir- 
acy between the distri bution of baryons and t he gravitational fiel d 
in spiral galaxies (e.g. jMcGaugh et alJuOOTI : iFamaev et alJuOOTT) . 
Indeed, it is baffling that such a simple prescription leads to uncan- 
nily successful pr edictions for galaxies ranging over five deca des 
in mass (see, e.g., ISanders & McGaughll200a Bekensteinll2006l for 
reviews), including our own Milky Way JFamaev & BinnevluOOa : 
iFamaev et aT]|2007n . 

On the other hand, recent developments in the theory of 
gravity have also added plausibility to the case for MOND 
through the advent of Lorentz-covariant theories yielding this be- 
havior in the ultra-weak gravity limit by means o f a dynam- 
ical (not necessarily) normalized f our-vector field (Bekenstein! 
2004IZlosnik et alj2007l:lBruneton & Esposito-Farese 2007; Zhaoj 



ABSTRACT 

Strong gravitational lensing by galaxies in Modified Newtonian Dynamics (MOND) has until 
now been restricted to spherically symmetric models. These models were able to account 
for the size of the Einstein ring of observed lenses, but were unable to account for double- 
imaged systems with collinear images, as well as four-image lenses. Non-spherical models 
are generally cumbersome to compute numerically in MOND, but we present here a class of 
analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged 
systems. We use them to obtain a reasonable MOND fit to ten double-imaged systems, as well 
as to the quadruple-imaged system Q2237 -1-030 which is an isolated bulge-disc lens producing 
an Einstein cross. However, we also find five double-imaged systems and three quadruple- 
imaged systems for which no reasonable MOND fit can be obtained with our models. We 
argue that this is mostly due to the intrinsic limitation of the analytic models, even though 
the presence of small amounts of additional dark mass on galaxy scales in MOND is also 
plausible. 
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and still being a far cry from a fundamental theory underpin- 
ning the MOND paradigm, these new theories r emarkably al 
low for new predictions regarding cosmology (e.g., ISkordis et al. 
200] 

2007^ and gravitational lensing (e.g 
!^95,; Zhaoetal. 2006; Angus et al 



new predictions regaraing cosmology (e.g., aKorais et al 
odelson & Liguorill2006l ; ISchmidt et alj|2007l; IZhang et al 



Chiu et al.ll2006l; lOin et al 



20071; IFamaev et al 



20071; 



Takahashi & Chiball2007l ; lFeix et alj2008 

In these Lorentz-covariant theories, gravitational lensing has 
actually been shown to work exactly in the s ame manner as i n 
General Relativity (GR) for static systems (e.g.. lBeken stein 2004J), 
except for the relation between the de nsity and the grav itational 
potential. This allowed, for instance, IZhao et alj ( 120061) to test 
MOND against the CfA-Arizona Space Telescope Lens Survey 
(CAST LES) data of image-splitti ng lenses, building o n the earlier 
work of loin et al.l (1995; also see lFerreras et al.ll2008r) . They were 
nevertheless restricted to models of spherical geometry, and were 
thus only able to account for the size of the Einstein ring of ob- 
served lenses, but not for the exact position of collinear images in 
double-imaged systems, and of course not for quadruple-imaged 
systems. This intrinsic limitation is due to the fact that the MON- 
Di an acceleration gM is related to the Newtonian one, ^^ = -VOiv> 
by JBekenstein & Milgron]|l984l) 



20071 ; ISanderil2005l ; [SkordisI |2008|) . Although rather fine-tuned 



ao / 



(1) 
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where 5 is a solenoidal vector field determined by the condition 
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that gn can be expressed as the gradient of a scalar potential O^. 
The function /j, controlling the modification of Newton's law, has 
the following asymptotic behavior: 



/i(x) ~ X X « 1, 

//(x) ~ 1 X » 1. 



(2) 



The field 5 is typically non-zero in non-spherical geometry, and 
calculating it usually requires the use of a nu merical solver for 
the non-li near field e quation of MON D (e.g.. ICiotti et alj|2006l ; 
iTiret & Combes 20071 : iFeix et alj2od3) . 

In this paper, we will demonstrate how to create simple ana- 
lytic models of non-spherical lenses with 5 = in MOND, choos- 
ing the covariant framework of tensor- vector-scalar gravity (TeVeS, 
lBekensteial2004l) . Without resorting to a numerical Poisson solver, 
these analytic models can thus be used to fit image positions in 
double-imaged and quadrupled-image systems. We present our an- 
alytic models in Sect. 2, devise a fitting procedure in Sect. 3, and 
finally expose and discuss the results in Sects. 4 and 5. 



2 THE HERNQUIST-KUZMIN MODEL 

2.1 Potential-density pair 

The Kuzmin disc JKuzminll 19560 . defined by a Newtonian gravita- 
tional potential of the form 

-GM 



<^N.K 



Vx2 + >- + (Izl + bf- 



h>0 



(3) 



is a well-known non-spherical case which satisfies S = in Eq. (TJ: 
For z > 0, Eq. ^ corresponds to the Newtonian potential generated 
by a point mass located at (0, 0, -b), in case of z < it turns into 
the Newtonian potential of a point mass located at (0, 0, b). Thus, 
above and below the disc, we eff^ectively have a spherical Newto- 
nian potential, which implies that truly 5 = in Eq. (TJ. 

Hereafter, the idea is simply to model lens galaxies by replac- 
ing the auxiliary point lens potenti al of the Kuzmi n disc with an 
auxiliary Hernquist lens potential ( lHernquistlll99Cl) : we shall re- 
fer to this model as the Hernquist-Kuzmin (HK) model. A simi- 
lar approach, using Plummer's model and a smooth transition at 
z = instead, leads to th e Plummer-Kuzmin model derived by 
iMivamoto & Nagail l ll975l) which provides a qualitatively good fit 
to the mass profile of observed galaxies. Although our proposed 
model is not a very good description of real galaxies, it enables us 
to derive fully analytic lens models in the context of MOND (see 
Sect. 12. 21 ) and to study the influence of non-sphericity on the ability 
to fit image positions. 

The Newtonian potential of the Hernquist-Kuzmin model 
takes the form 



9. 



-GM 



Vx' + y2 + (Izl + b)" + h 



(4) 



with b being the Kuzmin parameter and h denoting the core ra- 
dius of the Hernquist profile. Choosing different ratios h/b, this 
model will produce different Hubble type galaxies, going from a 
pure Kuzmin disc galaxy for h/b — > to a pure Hernquist sphere 
for h/b — > oo. To clarify this situation and to characterize the non- 
sphericity of the model, one may simply expand the r.h.s of Eq. l[4j 
far away from the disc (r^ = x^ +y^ + z')'- 



9^ 



-GM 

r + h 



\z\b 
(r + h)r 



-0{b^) 



(5) 



Using Poisson's equation, we find that the underlying density 
distribution is given by {R- = x^ + y-) 



Mh 



Phk 



In ^E? + (Izl + bf- ( V^' + (Izl + by + hf 



(6) 



The corresponding density contours in the (R, z) plane are plotted 
in Fig.[T]for different values of h/b. 

Considering the Hernquist-Kuzmin model for gravitational 
lensing, we choose the z-axis such that it is parallel to the line- 
of-sight and x,y are the Cartesian coordinates spanning the lens 
plane. Because we need to account for different possible orienta- 
tions of galaxies, we additionally have to rotate the disc. Defining 
(r'f = (x'f + iy'f + (Iz'l -1- bY, where 



X = (x cos - y sin (/>) cos 6* - z sin 9, 

y = xsinip + y cos 0, 

z = (x cos - y sin 0) sin 6 -I- z cos 0, 



(7) 



the angle (;r/2) - 6 being the inclination of the galaxy's symmetry 
plane w.r.t. the line of sight and (p the galaxy's position angle (PA), 
Eq. l|4) turns into 



*« 



-GM 

r + h' 



(8) 



2.2 Lensing Properties 

In this work, we consider MOND within the framework of TeVeS 
CBekenstein 2004), henceforth using units with c = 1 . Since most of 
the light bending still occurs within a small range around the lens 
compared to the distances between lens and source and observer 
and source, gravitational lensing works exactly the same way as 
in GR if the Newtonian potential is substituted by the total non- 

relativistic modified pote ntial. 

Then, according to ISchneider et al.l ([1993), the resulting de- 
flection angle can be written as 



2 fvxOdz 



(9) 



where (b is the total gravitational potential, Vx denotes the two- 
dimensional gradient operator perpendicular to light propagation 
and integration is performed along the unperturbed light path 
(Bom's approximation). Positions in the source and the lens plane 
are related through the lens equation: 



rj- 



-^f-A.,c?(f). 



(10) 



Here f = (x,y), ff is the 2-dimensional position vector in the source 
plane, and D^, Di, and D/j are the (angular diameter) distances be- 
tween source and observer, lens and observer, and lens and source, 
respectively. Furthermore, it is convenient to introduce the deflec- 
tion potential T(0): 



^¥0) = 2D 



f^iDi^, 



z)dz. 



(11) 



where we have used 9 = i/Di and D = DiJ(DsDi). If a source 
is much smaller than the angular scale on which the lens proper- 
ties change, the lens mapping can locally be linearized. Thus, the 
distortion of an image can be described by the Jacobian matrix 



Me) 



d/3 



■ K-Jl 

-72 



-72 

\- K + yi 



(12) 
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Figure 1. Contours of equal density in tlie (R,z) plane for the HK model 15) when h/b = 0.1 (top right), h/b = 1 (bottom left) and h/b = 10 (bottom right). 
Contour levels are (0.01, 0.003, 0.001, ...)M/b^ (top right); (0.001, ...)M/b^ (bottom left); (0.0003, ...)M/b^ (bottom right). The top left panel illustrates the HK 
model: At the point {R, ~\z\) below the disc, the potential Eq. (4) is identical with that of a Hernquist distribution whose origin is located at a distance b above 
the disc's centre. 



with the angular coordinate /3=fj/Ds. The convergence k is directly 
given by the deflection potential T: 



2 " 



(13) 



Similarly, the shear components yi, 72 can be calculated from 



1 /a^'P 



71 = 5 



5(9; 36: 



a-f 

^ 



72 



36,302 



(14) 



r 



^4 



Ti+n- 



Due to Liouville's theorem, gravitational lensing preserves the sur- 
face brightness, but it changes the apparent solid angle of a source. 
The resulting flux ratio between image and source can be expressed 
in terms of the amplification A, 



(l-Kf-y- 



(15) 



Consequently, the flux ratio between two images A and B is f/^g = 
A^/Ag. Points in the lens plane where A"' = form closed curves, 
the critical curves. Their corresponding curves located in the source 
plane are called caustics. Images near critical curves can signifi- 
cantly be magnified and distorted, which, for instance, is indicated 
by the giant luminous arcs formed from source galaxies near caus- 
tics. 

Due to the deflection by the gravitational potential, light rays 
traveling from a source to an observer at redshift z = will be 



delayed in time: 



m = 



l+z, 



^(l(6-py--^'{e)Y (16) 

where zi is the lens' redshift. If the deflection potential is known, 
Eq. l |16t allows to calculate the relative time delay between differ- 
ent images. 

Choosing a certain smooth transition from Newtonian dynam- 
ics to MOND (see Zhao et al. 2006, Sect. 8.2 for details) and as- 
suming spherical symmetry, we have the following relation for the 
total acceleration in TeVeS: 



giuir) = gN(r) + ^jgN{r)aQ, 



(17) 



where Aq = 1.2 x lO'^ms"^. Exploiting the above and introducing 
Zo = (xcos(/i - 3; sin 0) tan (?, the deflection angle's x-component 
yields 



a J = 2(x - b cos (f) cos 6) 



+ 2{x + b cos cos (?) 






GM 

{/ + hy 

GM 

(r + hf 



^/GM^ 



ao 



r +h 



^JGMao 
r + h 



(18) 



The integral l llSt can be evaluated by means of elementary calculus, 
but as the resulting expression is quite lengthy, we shall skip its 
presentation at this point. Analogously, the closed analytic form 
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for ay can be derived, and as as a consequence, this is also true for 
the lensing quantities k and 7. 



2.3 Cosmology and distances 

Concerning the calculation of distances in gravitational lensing, 
we shall adopt a standard flat ACDM cosmology with CI,,, = 0.3 
and SIa = 0.7. This choice is justified by the fact that many co- 
variant formulations of MOND mimic the behavior of a ACDM 
model, accounting for marginal differences tha t will have no si 
nificant impact on our analysis. For instance, IZhao et al.l j200i 
presented a minimal-matter open model with a cosmological con- 
stant whi ch works to good accuracy for redshifts ; < 3. Fur- 
thermore, lAngus et all ( 120071) proposed the flat yuHDM (MOND 
Hot Dark Matter) cosmological model based on the assumption of 
massive (sterile) neutrinos. If provided with a cosmological con- 
stant, this model is quite similar to the standard flat ACDM model, 
being capable of e xplaining the cosmic microwave background 
JSkordis et al.l200q) . Recent work has also shown that it seems pos- 
sible to construct certain covariant formulations of MOND gi ving a 
cosm ology identical to ACDM in the matter-dominated era jZhad 
l2007h . 



tions with an anomalous mass-to-light ratio or a large flux anomaly: 

h 



P 



(log FP}- + (cos if + 
M 



'r 



l°g xobs 



b + h 

2 



log 



M,, 



(21) 



The deviation from the Fundamental Plane is measured by 
log FP = log(h/hi ) - l.26log(M/Mi) and hi = 0.72kpc and 
M, = 1.5 X 10"Mo JChen & Zhaoll2006l) . 

We choose a very small regularization parameter A ~ 
(0.003")^, and minimize the following regularized "^--like" quan- 
tity. 



77^ = A? + A^ + AP, 



(22) 



for 14 double-imaged systems and four quadruple-imaged systems 
in CASTLES. Note that we also check that our results are insensi- 
tive to the detailed choice of the regularization parameter and that, 
due to the sufficient amount of constraints (position of lens and 
images), the fitting procedure for quadruple-image lenses is per- 
formed with A = 0. The results are shown in Table [T] and Table |2] 
respectively. Finally, note that the ob served mas s of ea ch lens was 
calculated according to Sect. 7.1 of IZhao et al.l ( 120061 Eqs. (73)- 
(75)). 



3 FITTING PROCEDURE 



We will follow iKavseii jl990l) to model the lens systems: For each 
pair of images / and j, when tracing one light-ray back for each 
observed image to the source plane, the source position obtained 
from Eq. dlOl l should be the same for both images. We can thus 
simply compare the resulting source position for each image by 
computing their squared deviation. 



Aj = 2 {(x,i - x,jf + (y„ - y,jf) , 



(19) 



where x, and Vj is the source position in Eq. ((TO}. This is a mea- 
surement of how well the images retrace back to a single point in 
the source plane. 

Another quantity to minimize is the deviation of the lens cen- 
tre from the observed optical centre, given by 



A2 = [(x,f + (y,)-) , 



(20) 



However, our model has generally nine fitting parameters (the 
lens mass M, the Kuzmin length b, the Hernquist length h, the PA 
angle, the inclination 7, the source position Xj,yj, and the lens posi- 
tion x/,y/), while for a double-imaged system we have only four 
constraints from the two image positions, and another two con- 
straints from the observed lens optical centre. The problem is thus 
ill-posed. 

To cure this, and ensure the uniqueness of the solution, we use 
a regularization method by adding a regularization term in the min- 
imization. This term is penalizing solutions deviating from the Fun- 
damental Plane as well as face-or|j and discy solutions, and solu- 



4 FITTING RESULTS 

4.1 Double-imaged systems 

Setting A, < 0.01" as a reasonable threshold for acceptable fits of 
the HK lens. Table [T] shows that our model is able to describe the 
observed image positions of all double-imaged systems, with quite 
a number of these systems yielding plausible parameters within the 
context of MOND/TeVeS. Additionally, the HK model seems to be 
able to explain the flux ratios of these binaries in almost every case. 
However, there are a few outliers which we will discuss in the 
following. Since the model should be capable of reproducing all ob- 
servational constraints and the lens mass should have a value close 
to the stellar mass (M/M, ^ 1) in MOND/TeVeS, these are charac- 
terized by very poor fitting parameters in terms of large differences 
between predicted and observed flux ratios and/or anomalous mass 
ratios M/M, (deviation larger than a factor of 3). 



4.1.1 RXJ0921+4529 

The system RXJ0921+4529 contains two Zs = 1.66 quasars and 
a H = 18.2 spiral galaxy l ocated in between the quasar images. 
According to lMunoz et alj j200lh . this galaxy lens is quite likely 
to be a member of a. zi = 0.32 X-ray cluster centreed on the ob- 
served field. Clearly, RXJ092 1-1-4529 does not conespond to an 
isolated system, which complicates the situation in TeVeS and pro- 
vides a possible explanation for the extremely poor fit/mass ratio 
(M/M, » 59). The presence of a cluster could have caused diffi- 
culties in fitting the lens as the impact of non-linear effects may 



' As there is strong observational evidence supporting that the sys- 
tem B0218+357 corresponds to a nearly face-on spiral galaxy (e.g. 



iGrundahl & Hiorth|[T99a lYork et al.ll2005l) . we choose the regularization 
term for this particular lens such that edge-on solutions are penalized in- 
stead. Further relaxing the penalties w.rt. the Fundamental Plane and the 
observed flux ratio in Eq. 1211 . the fit substantially improves, correspond- 
ing to a factor of 20 in A,. 



^ In case of RXI092 1+4529, however, our choice of A creates an over- 
regularization effect, which results in a best-fit lens ma ss that is roughly by 
a factor 10 smaller than estimated bv lzhao et alj 1200g ) fitting the system's 
Einstein ring size. Decreasing the regulaiization parameter to /I ~ (3 X 
[0-4// j2 jj, gjjjg jQ [-esQive this issue, with A, dropping by a facto r 10 and the 
lens m ass now being in accordance with the previous estimate of lzhao et alj 
i200d) (see TableQ. 
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Table 1. Fitting results fo r selected two-image lens systems from the CASTLES sample: In the table, the observed lens mass M, is calculated accord ing to Sect. 
7.1 oflzhao et alj 120061 Eqs. (73)-(75)), the parameter r;, is the Hernquist length expected from half-light measurements (values are taken from Izhao et alj 
120061) . We do not give tj^, but instead we list A, and compare the inferred values of PA, the Hernquist length h, mass and flux ratio to observations. Additionally, 
we predict inclination and time delay for the particular lens models. Outliers are characterized by large differences between predicted and observed flux ratios 
and/or anomalous mass ratios M/M, (deviation lai'ger than a factor of 3) like, for instance, in case of RXJ092 1+4529 which resides in a cluster. Note that the 
fitted lens position is given by (xi,yi) » (0, 0) for all lenses. 



Lens 


Zl 


h/h 


h/n, 


M(fit/obs) 


PA 


inch. 


A 




/As(fil/0bs) 


St(fit/obs)) 








kpc 


10" Mo 


degree 


degree 


arcsec 




days 


QO 142- 100 


0.49 


0.25 


1.34/1.6 


1.70/4.08 


72.2 


90.0 


2.43 X 


10-* 


8.06/8.22 


151.5/- 


B02 18+357 


0.68 


1.0 


2.19/1.8 


2.69/2.67 


-22.6 


6.94 


7.75 X 


10-^ 


0.759/0.587 


7.52/10.5 


HE0512-3329 


0.93 


0.24 


1.45/1.8 


1.49/2.9 If. 


28.2 


90.0 


3.90 X 


10-" 


0.0013/1.175 


19.6/- 


SDSS0903+5028 


0.39 


0.76 


1.83/1.8 


2.77/3.80 


-30.4 


90.0 


9.90 X 


10--* 


2.29/2.17 


135.2/- 


RXJ092 1+4529 


0.31 


0.037 


7.59/1.8 


20.0/0.34 


60.2 


90.0 


5.85 X 


10-* 


3.623/3.591 


167.2/- 


FBQ095 1+2635 


0.24 


0.13 


1.20/0.32 


0.47/0.31 


60.3 


90.0 


1.23 X 


10-* 


2.74/3.53 


13.2/- 


BRI0952-0115 


0.41 


0.055 


2.20/0.29 


0.58/0.27 


124.1 


90.0 


5.04 X 


10-* 


3.52/3.52 


8.11/- 


Q0957+561 


0.36 


1.55 


1.21/5.23 


6.94/8.44 


40.0 


90.0 


1.97 X 


10-^ 


14.3/1.08 


752.4/417.0 


Q1017-207 


0.78 


0.0092 


2.39/1.19 


0.83/0.74 


88.8 


89.9 


2.16X 


10-* 


0.73/0.72 


29.0/- 


B 1030+071 


0.60 


0.10 


0.84/1.50 


1.85/1.66 


29.3 


90.0 


8.04 X 


10-^ 


36.6/36.6 


346.8/- 


HE 1104-1805 


0.73 


0.33 


0.58/2.48 


4.91/3.32 


61.9 


90.0 


1.96 X 


10-^ 


0.35/3.85 


321.2/- 


B 1600+434 


0.41 


0.18 


1.64/1.8 


1.01/0.40 


36.8 


90.0 


2.09 X 


10-* 


0.83/0.84 


32.2/51.0 


PKS 1830-2 11 


0.89 


0.48 


2.75/1.8 


1.33/1.48 


62.3 


90.0 


4.14 X 


10-* 


157.3/157.3 


32.7/26.0 


HE2 149-2745 


0.50 


0.026 


0.94/11.4 


1.04/2.00 


-30.0 


90.0 


2.30 X 


10-* 


6.53/4.19 


90.7/103.0 


SBS0909+523 


0.83 


0.19 


3.02/1.8 


2.92/13.5^ 


49.2 


90.0 


1.84 X 


10-^ 


1.42/1.42 


65.9/- 


hatlZhao et alj i2006h used 


i different value for M, 


based on a wrong magnitude in an 


older version of the CASTLES data set. 



be important. In addition, note that there are still s ome unresolved 
issue s in MOND/TeVeS concerning clusters ( e.g. ISanderslll999l 
I200I I200I lAngus et al.ll2007l ; iFeix et al.ll2008l) . 



4.1.2 Q0957+561 

The gravitational lens Q0957+561 is the most thoroughly studied 
one in literature. The system involves a radio-loud quasar at redshift 
Zs = 1.41 which is mapped into two images by a brightest clus- 
ter galaxy (BCG) and its parent cluster at redshift zi = 0.36 (e.g. 
iKeeton et al.ll200(]| : lBarkana et al.ll 19991 and references therein). It 
is also known that the lens galaxy has a small ellipticity gradient 
and isophote twist which are properties the simple HK model can- 
not account for. Together with the fact that the lens is embedded 
into a cluster, this might be a reason for the huge discrepancy be- 
tween observed and predicted flux ratio in the context of modified 
gravity. 



4.1.3 HE1104-1805 

The lens galaxy's colors are in agreement with a high-redshift 
early-type galaxy, and its redshift is roughly estimated as zi = 0.77 
JLehar et ahluOOOl and references therein). Concerning its lensing 
properties, the system HE 1 104-1805 is quite uncommon in a sense 
that the lens is closer the bright image, rather than the faint one. 
As is known from lensing within the standard GR/CDM paradigm, 
simple models can create such configurations only for a narrow 
range of parameters due to the peculiar flux ratio. Assuming sim- 
ple ellipsoidal lens models, however, these parameters imply a large 
misalignment between the light and the projected density. The only 
possibility to align the mass with the light is to have a shear field 
being approximately twice as strong as estimated from the particu- 
lar lens model. 

Furthermore, the observed image separation is by a factor 2-3 



larger than that of a typical lens, strongly suggesting that the sep- 
aration is enhanced by the presence of a group or a cluster So far, 
however, there has been no direct observational evidence for such 
a structure in the lens' surrounding area. 

Analog to the above-mentioned lens systems, the unsatisfying 
fit and the correspondingly inferred flux ratio might be a result of 
both lens environment and model limitations (see also Sect. l4.3l and 
W3. 



4.1.4 SBS0909+523 

SBS0909+532 shows two images of background quasa r source at 
Zs = 1.377 separated by 1.11" (e.g. iLehar et alJuOOOl and refer- 
ences therein). Optical and infrared HST images indicate that the 
lensing galaxy has a large effective radius and a correspondingly 
low surface brightness. Additionally, the lens galaxy's redshift is 
estimated as zi = 0.83 JLubin et al.ll2000h . and its total magnitude 
in the //-band has been measured as H = 16.75 ± 0.74. Although 
the lens galaxy's colors are poorly measured, they seem consistent 
with those of an early-type galaxy at the observed redshift. 

The large uncertainties are a r esult of the difficu lty in subtract- 
ing the close pair of quasar images JLehar et al.l2000l) . For instance, 
the uncertainty in the /-band magnitude, / = 18.85 ± 0.45, allows 
for a deviation of the mass estimate M, by a factor of rou ghly 2.3 on 
a 2cr l evel, where we have again used Eqs. (73)-(74) of lZhao et al.l 
1 I2OO6I) . Thus, we argue that the low mass ratio (listed in Table [TJ 
may be entirely due to these uncertainties in observed magnitudes, 
with better constrained observations possibly softening the problem 
in MOND/TeVeS. 



4.1.5 HE0512-3329 

The system HE05 12-3329 was discovered as a gravitational lens 
candidate in the course of a snapshot survey with the Space Tele- 
scope Imaging Spectrograph (STIS), with the images of the quasar 
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Figure 2. Critical curves (black) and caustics (red) of the best-fit Hemquist- 
Kuzmin model for Q2237+030: The empty and filled squares denote the 
observed positions of images and source, respectively. 
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Figure 3. Convergence map k of the best-fit Hernquist-Kuzmin model for 
Q2237-h030, with the outer contour level being at k = 0.7 and increasing in 
steps of 0. 1 up to a level of *: = 2.0. 



source being separated by 0.644" jWucknitz et al.ll2003l . and ref- 
erences therein). Although the lensing galaxy has not been de- 
tected yet, measurements of strong metal absorption lines at red- 
shift z = 0.93, identified in the integrated spectrum, hint towards a 
damped Lyo- system intervening at this redshift. 

Analyzing sepa rate spectra of both image components, 
IWucknitz et al.l j2003h point out that both differential extinction and 
microlensing effects significantly contribute to the spectral differ- 
ences and that one cannot be analyzed without taking into account 
the other. For lens modeling purposes, the observed flux ratio can 
therefore only be used after correcting for both effects. Thus, the 
large discrepancy between predicted and observed flux ratio might 
be a consequence of neglecting the above mentioned effects, rather 
than being intrinsic to MOND/TeVeS. 



4.2 Quadruple-imaged systems 

As we can see from Tabled most of the quadruple-imaged systems 
are very poorly fitted by our analytic HK lens model for MOND. In 
accordance with our goodness-of-fit criterion (A, < 0.01") intro- 
duced in Sect. 14. II there is just one system where the model is able 
to predict the image positions in a satisfying manner. Additionally, 
none of the observed flux ratios can be explained. 

The only acceptable fit is given for Q22 37-I-030, the nearby 
Einstein cross (zi = 0.04; iHuchra et alj|l985l) . which is the only 
true bulge-disc system in our set. Also, its physical Einstein ring 
size in the lens plane is very small, R^ « 0.7kpc (in B 14 22-1-231, for 
instance, it is already by a factor of roughly 10 larger: iRusin et al.l 
120031) . Nevertheless, it is not possible to give a reasonable explana- 
tion for the flux ratios using our smooth MOND lens model. Taking 
effects due to microlensing into account, which are not considered 
in this work, could be able to relax the situation. Note that the lens 
galaxy actually contains a bar feature JYedl 19881) which is ignored 
in our analysis. 



4.2.1 PG1115+080 

The lens galaxy in PGl 1 15-1-080 and its three neighbors belong to 
a single group at z; = 0.31 1, with the group being centreed south- 
west of the lens galaxy's position l IKundic et al.lll997lMlmDev et al.l 
119981 and references therein). Reasonable fits of this lens typically 



involve a significant amount of external shear in the context of 
GR/CDM. Moreover, the observed anomaly of the flux ratio (~ 0.9) 
between two of the images strongly hints towards an additional per- 
turbation of the system caused by a satellite galaxy or a globular 
cluster Similar to Sect. 14.11 we have a gravitationally bound sys- 
tem which will involve a more complicated approach in TeVeS than 
provided by the isolated HK model. 



4.2.2 B1422+231 

The system B 1422-1-231 shows al most the same characteristics as 
PGl 115-1-080 JKundicetalll 19971 , and references therein). Again, 
the lens belongs to a gal axy group which is cent reed south of the 
lens galaxy (zi = 3.62). iKormann et alj ( Il994ah were able to fit 
the l ensing system using a very flat singular isotherm al ellipsoid 
(SIE: lKassiola & Kovneiill993l : iKormann et alj|l994bl) plus an ex- 
ternal shear field. However, HST observations revealed that the lens 
galaxy's optical axis ratio is much closer to unity than assumed 
for the flat SIE, favoring rounder lens models with larger external 
shear. 



4.2.3 SDSS0924+0219 

Estimated colors and magnitudes of the lens galaxy ar e consistent 
with those of a typical elliptical galaxy at zi = 0.4 jlnada et al.l 
I2OO3I and references therein). Although the lens environment does 
not show any nearby objects perturbing the system, quite an amount 
of external shear is needed to obtain a satisfying fit to observa- 
tions, with the lens being typically modeled by a (flattened) singular 
isothermal sphere (SIS). Additionally, microlensing plays an im- 
portant role in explaining the observed flux ratios within GR/CDM, 
which is likely to be true in TeVeS as well. 



4.3 Maximum non-spherical shear of a Kuzmin lens 

As we have seen, the outliers in our selection of quadruple- 
image lenses correspond to systems with a large external shear. 
In PG1115-f080, for example, this is due to a neighboring galaxy 
group. However, the same situation also appears in uncrowded en- 
vironments, usually constraining the lensing potential to require 
a substantial ellipticity. From Sect. 14.21 it seems that our present 
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Table 2. Fitting results for selected four-image lens systems from the CASTLES sample: Note that all positions (RA and Dec.) are given in units of arcsec. The 
observed PA and inclination of Q2237+030 (major-axis) are PA= 77.2° and (' = 64.5°, respectively, assuming a circular face-on disc. Replacing the auxiliary 
Hernquist with a Jafl'e profile barely changes the numbers: inclination and PA change by about 5 degrees, the predicted mass by roughly 10 per cent. 





PGl 115+080 


Q2237+030 


B1422+231 


SDSS0924+0219 


Zl 


0.31 


0.04 


0.34 


0.39 


Zs 


1.72 


1.69 


3.62 


1.52 


D, (kpc) 


957.2 


163.6 


1020.2 


1116.6 


D, (kpc) 


1874.2 


1874.0 


1637.6 


1867.0 


D,, (kpc) 


1413.2 


1810.8 


1341.7 


1252.1 


Image A 


(-0.947, -0.690) ± 0.003 


(-0.075, -0.939) ± 0.003 


(0.375, 0.973) ±0.003 


(-0.162, 0.847) ±0.003 


Image B 


(-1.096, -0.232) ±0.003 


(0.598, 0.758) ±0.003 


(0.760, 0.656) ±0.003 


(-0.213, -0.944) ±0.003 


Image C 


(0.722, -0.617) ± 0.003 


(-0.710, 0.271) ±0.003 


(1.097, -0.095) ±0.003 


(0.823,0.182) ±0.003 


Image D 


(0.381, 1.344) ±0.003 


(0.791, -0.411) ±0.003 


(-1.087, -0.047) ±0.003 


(-0.701, 0.388) ±0.003 


Source 


(-0.011,0.091) 


(0.027,-0.0051) 


(0.089, 0.030) 


(-0.024, -0.047) 


Lens 


(-0.0011,-0.0041) 


(0.00066,0.00096) 


(-0.00093, 0.0065) 


(0.019,-0.0051) 


M(Mfi,/M.)(xlO"Mo) 


7.80/1.23 


0.78/1.19 


4.83/0.77 


2.80/0.32 


h (kpc) 


2.25 


0.44 


8.42 


1.57 


b/h 


0.56 


1.85 


0.29 


2.17 


PA angle (degree) 


244.8 


246.6 


117.9 


266.4 


Inclination (degree) 


44.5 


30.6 


48.6 


40.5 


Aj (arcsec) 


0.0402 


0.0026 


0.0593 


0.0612 


Flux ratio(obs) 


4.03:2.53:0.65:1 


2.62:1.64:1.30:1 


31.1:34.6:18.4:1 


12.5:5.68:4.81:1 


Flux ratio(fit) 


3.98:4.15:1.40:1 


0.81:0.66:0.68:1 


8.56:6.53:7.51:1 


1.66:0.69:0.86:1 



analytic model is not able to generate such a potential in most 
cases. As is known, almost all quadruple-image d systems show ev- 
idence for the need of an external shear fie ld jWitt & MadI 19971 ; 
iKassiola & Kovneiil993l : lKeeton et all 19971) by violating a certain 
inequality of the image positions. It is perhaps not surprising that 
the current isolated HK model fails to fit these lense^ 

To gain a better understanding about this issue, we consider a 
pure edge-on Kuzmin lens {h = 0) and derive the maximum varia- 
tion of the shear at the Einstein radius Rg by comparing its values 
on the major and minor axis. For this reason, let us introduce a 
quantity Q being given as follows: 



Q 



y(RE,O)-y{0,RE) 
y{RE,0) + j(0,REy 



(23) 



The parameter defined above will indicate the level of the shear 
field's non-sphericity at the Einstein radius and is a function of the 
dimensionless radius RE/b. Note that, in case of the Kuzmin lens, 
the quantity Q depends on redshift. 

Figure |4] shows Q as function of RE/b for the pure Kuzmin 
model (solid line), assuming aoD/c^ = 0.03. This value has been 
chosen in accordance with the majority of lens redshifts in the 
CASTLES sample, and changing it does have no significant qual- 
itative impact on the basic outcome. Additionally, we also present 
the result for the SIE model l IKormann et al.l 19940) . with the poten- 
tial axis ratio varying from 0.7 to 0.9 (shown by horizontal lines). 
As we can see, the Kuzmin model becomes comparable to a very 
round SIE if Re /b> 10. 

To obtain a sufficiently strong quadruple moment/non- 
spherical shear at the Einstein radius (Q > 0.2), these disc-only 
models must satisfy the condition b > 0.2Re. In case of PGl 115, 
the observed ring size can be estimated as Re ~ 5kpc, so to fit 




RE/b 

Figure 4. Non-spherical shear parameter Q for a simple TeVeS Kuzmin lens 
(solid line), assuming oqD/c^ = 0.03. Additionally, we show the results for 
a SIE model with a potential axis ratio of 0.9 (dot-dashed line), 0.8 (dotted 
line) and 0.7 (dashed line), respectively. 



four images, one might actually expect that b > r x Ikpc. How- 
ever, trying to fit the above mentioned Einstein ring size, using the 
stellar mass only, we also find that this would need a Kuzmin pa- 
rameter close to zero {b a 0), corresponding to a very concentrated 
point-like lens. Although we have only given a plausibility argu- 
ment, rather than a rigorous proof, this could explain why we can- 
not find a value of b that meets both requirements and why the HK 
model mostly fails to fit quadruple-imaged systems. 



^ Note that our analysis does not take into account external shear effects, 
which would complicate the relation between lens mapping and associated 
density distribution due to non-linearity in modified gravity. While the main 
task of this paper is to explore the capability of the HK model, such contri- 
butions should certainly be addressed in future work. 



4.4 Experimenting with hypothetical lenses 

Another possibility of investigating the fitting capability of our 
model is to generally explore its parameter space and to study the 
structure of critical curves and caustics. To avoid any limitations 



8 HuanYuan Sharif Martin Feix, Benoit Famaey and HongSheng Zhao 



a=Q. b/(h+b)=6.36-, ,incli=44.5 



' ' 


1 


1 1 1 _ 1 


.11 


- 




„= 


-1, b/(h-b)=0,36, in 


=li=44,5- 












' 


- 








- 


_ 


" 


- . _ ■ ' ° 




_ 




^i* 1. i 






: 







1 1 1 1 




1 , 


- 


1 




a = l, b/(h+b) 


=0.36. 


ncli = 60 


- 




1 


- 


l_l L 1_, 1 









-2 1 Q, 1 2 -2 -1 Q, 1 3 -2 -1 Q, 1 2 

Figure 5. Experimenting with hypothetical lenses: Shown are the critical curves (black) and caustics (red) of different Dehnen-Kuzmin models characterized 
by the parameters a, h/{h + h) and inclination ("core" radius and PA are fixed to h = 0.72kpc and 11.2° , respectively). All models assume a lens mass of 
M = 8 X 10"Mo which is approximately 8 times the stellar mass of the lens galaxy in PGl 1 15+080. The empty and tilled squares denote the observed image 
and source positions of PGl 1 15+080. 



o=l, b/(h+b)=0,3S, incli=44.5 



a = l, b/(h+b) = 0.36, incli=44.5 
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b=1.14 arcsec, 7 = 07. 8 =B8.7 
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Figure 6. Effects of reducing the lens mass: Shown are the critical curves (black) and caustics (red) of a Dehnen-Kuzmin model (a = 1, bl(h + h) = 0.38, 
h = 0.72kpc, PA= 77.2° and; = 44.5°), assuming M = 8x 10"Mo (left panel) andM = 5x 10"Mo (middle panel), respectively. The empty and filled squares 
denote the observed image and source positions of PGl 115+080, the stellar mass of the lens is estimated as M, » 10"Mo. Right panel: Critical curves (black) 
and caustics (red) of the best-fit SIS+7„( model given by Eq. J25j. 



that might be due to the particularly chosen radial profile, we fur- 
thermore replace the aux iliary Hernquist profile with the more gen- 
eral model proposed by lOehnenl jl993h . Its Newtonian potential 
and the corresponding density profile read as 



<^N.D(r) ■ 



Pair) ■■ 



GM 



-1 + 



[r + hl \' 



Mh(l + a) 



(24) 



where h is a characteristic length of the model. Depending on the 
value of a, the Dehnen model represents different density distri- 
butions, ranging from quite cuspy to more broadened profile s. For 
= and a = I, Eq. ( I24t reduces to the models of Jaffe Uaffd 
Il983h and Hernquist, respectively. 

Allowing different values for a, we repeat the fitting procedure 
for the quadruple-imaged systems discussed in Sect. |4.2| The result 
is basically the same as for the HK model, with the parameters 
listed in Table [2] not significantly changing. In case of the Jaffe 
profile (a = 0), for instance, inclination and PA are altered by about 
5° and the predicted mass by approximately 10%. 

To further illuminate this obvious insufficiency of our model, 
let us have a more detailed look at the caustic structure, taking the 
system PG1115+080 as an example: Choosing a plausible setting 
for the lens system in MOND, we fix its size to h = 0.72kpc and 



the PA to 77.2° (observed value). In accordance with the best-fit 
results, we additionally assume a lens mass of M = 8 x 10" Mq and 
vary the Dehnen index a, the model's "discyness" b/{h + b) as well 
as the inclination on a range from -I to 1, 0.1 to 0.9 and 10° to 
90°, respectively. For a selection of such lenses, the corresponding 
critical curves and caustics are shown in Fig. [S] Then, among all 
resulting lens models, we select those which exhibit the strongest 
(non-spherical) shear, corresponding to a large astroid caustic size. 
Since the lens mass should be close to the stellar mass (M/M, ^ 1) 
in MOND/TeVeS, the idea is now to stepwise decrease the mass of 
these models. In all cases, we find that, due to the caustics' con- 
traction, the source crosses the astroid caustic way before M/M, 
reaches unity, thus not corresponding to a quadruple-imaged sys- 
tem anymore. Typically, the crossing seems to take place when the 
lens model's mass is roughly around 4 - 6 x 1O"M0. For a = I, 
b/(h + b) = 0.38 and an inclination of 44.5°, this situation is illus- 
trated in the left and middle panel of Fig.|6] Note that we have kept 
the source position fixed at (-0.01 1, 0.091)" for our analysis, with 
the lens being centred at the origin. 

Again, this provides a possible explanation why the Dehnen- 
Kuzmin model (including the HK model) mostly fails to fit 
quadruple-imaged systems, supporting our earlier conclusion from 
Sect. 14. 31 Given that M/M, ^ 1 in MOND, our model is obviously 
not able to generate sufficiently strong shear (hence large caustics) 
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and a large Einstein ring at the same time. For comparison, we 
also present the resulting caustics and critical curves of a best-fit 
SIS+7jv( model in the right panel of Fig.|6] As is known, its deflec- 
tion potential can be expressed as 



'V(^,9) = c^+^fcos{2(e-9^)). 



(25) 



Choosing the lens' position (x,,y,) = (0.0028, 0.0048)", c = 1.14", 
7 = 0.07 and 0y = 88.7°, the above model is able to fit observations 
of PGl 115+080 satisfyingly. 



5 SUMMARY AND DISCUSSION 

In this paper, we presented a class of analytic non-spherical mod- 
els, called the Hemquist-Kuzmin models, that we applied to fit 
double-imaged and quadruple-imaged lens systems in the context 
of MOND, using the covariant framework of TeVeS. These ana- 
lytic models interpolate smoothly between a Hernquist sphere and a 
Kuzmin disc, and their great advantage is that the solenoidal field of 
MOND, which normally appears in Eq. ([TJ for non-spherical con- 
figurations, is zero. Note that these models are not, stricto sensu, 
bulge-disc models (see Fig.[T). 

After having worked out the lensing properties of the mod- 
els, we devised a fitting procedure in order to fit 15 double-imaged 
systems and four quadruple-imaged systems in the CfA-Arizona 
Space Telescope Lens Survey (CASTLES). Since the fitting prob- 
lem is ill-posed, especially in the case of double-imaged systems, 
we used a regularization method to ensure the uniqueness of the 
solution, by penalizing fits deviating from the Fundamental Plane 
as well as face-on and discy fits, and fits with an anomalous mass- 
to-light ratio or a large flux anomaly. 

From this fitting procedure, we find that our model is 
able to describe the observed image positions of all analyzed 
double-imaged systems, and we can safely conclude that 10 of 
these systems yield plausible parameters within the context of 
MOND/TeVeS. Additionally, our analytic model is able to explain 
the flux ratios of these binaries in almost every case. One can thus 
conclude that, in general, galaxy lenses do not need dark matter 
in MOND/TeVeS. Note that the implied masses for most of these 
lenses are quite si milar to tho s e deriv ed from the spherically sym- 
metric models of IZhao et al.l ( 120061) , but that the big advantage 
of our non-spherical model is its ability to fit the precise image- 
positions rather than just the size of the Einstein ring. 

On the other hand, 5 double-imaged systems do not provide a 
reasonable fit: while for two of these systems, the found problems 
are likely to be solved by considering observational uncertainties, 
a more accurate model and/or additional effects such as extinction 
and microlensing, the other three lenses appear to be lacking an 
obvious explanatiorO. It is however quite striking that all these re- 
maining outliers are actually residing in (or close to) groups or clus- 
ters of galaxies. This means that non-linear effects could have an 
important incidence, since lensing in MOND is much more sensi- 
tive to the three-dimensional d istribution of the lens (and of the 
environment) than in GR (e.g. iFeix et alj 120081) . Moreover, it is 
known for a while that additional dark matter is needed fo r clusters 
of galaxies in MOND (e.g. Sanders 1999. 2003. 2007 ; Angus et al.l 



* Note, however, that the stellar mass estimates depend on the adopted ini- 
tial mass function and stai' formation rate, and can vary by a factor of 4 in 
the R-band, which could partly solve the problem of the mass-ratio discrep- 
ancy, but not the flux ratio anomalies. 



l2007l : lMilgromll2008h , a nd it has recently b een shown that this was 
the case for groups, too JAngus et alj|2008l) . Possible explanations 
for this "cluster dark matter" in the context of MOND range from 
the presence of numerous clouds of cold gas JMilgrom 2008) to th e 
existence of sterile neutrinos with a 5eV mass JGentile et al.l2008l) . 
through the non-trivial effects of the vector field (or of an addi- 
tional scalar field) in covariant formulations of MOND other than 

I " ' ' " ^ 

TeVeS (e.g. [Sanders 2005; Zhao 2007"). In any case, many studies, 

including the recent analysis of the velocity dispersions of globular 

clusters in the halo of NGC 1399 JRichder et aLll2008l) . have also 

provided evidence for such dark matter on galaxy scales in MOND: 

this is typical for g alaxies residing at the centre of clusters only, and 

can be interpreted l IRichtler et al . 2008) as a small-scale variant of 

the aforementioned problem of MOND in clusters. This could thus 

provide an additional, physical, reason to the poor fits obtained here 

for the two-image lenses residing in groups or clusters. 

For the four quadruple-imaged systems, it is a different 
story: the only acceptable fit is obtained for the Einstein cross 
Q2237+030, but even in this case, the observed flux ratios can- 
not be reproduced. However, the anomalous flux ratios here are 
most likely due to microlensing effects which have not been con- 
sidered in our analysis. We can thus conclude that MOND does 
not provide a solution to the flux anomaly issue, mainly because 
smooth MOND models naturally predict smooth amplification pat- 
terns. Among the 3 very poorly fitted lenses, only PGl 1 15+080 and 
B 1422+231 appear in a crowded environment, which could cause 
the same perturbing effects as for non-isolated double-imaged sys- 
tems; the remaining lens, SDSS0924+0219, appears relatively iso- 
lated. 

We thus argue that, especially in this particular case, the poor 
fits are due to the intrinsic limitation of our analytic models: we 
indeed showed that our models were unable to produce a large 
Einstein ring and a large shear at the same time. Although we did 
not present a rigorous proof, we tried to make this limitation plau- 
sible by analyzing the maximum non-spherical shear of a TeVeS 
Kuzmin lens as well as the caustic structure of different Hemquist- 
Kuzmin models. To m ake the mass more concentrated, we also 
tried a Dehnen model ( lDehnenlll993n , but this did not provide a 
satisfactory solution either. Again, note that our analysis did not 
consider any contribution due to external shear effects. 

We therefore conclude that our analytic models provide good 
MOND fits to the image positions of isolated two-image lenses, 
but that some problems are encountered for non-isolated two-image 
lenses. On the other hand, we showed that our models were barely 
able to fit quadruple-imaged systems. The present study has thus 
pinpointed some lenses for which a full three-dimensional numeri- 
cal model should be devised in MOND. While these analytic mod- 
els did obviously not yet represent a definitive test of MOND with 
gravitational lensing, they nevertheless provided a new step in un- 
derstanding this quite unexplored research area, and in isolating the 
possibly challenging lensing systems for the future. 
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